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Abstract. Determination of the quasi-elastic scattering cross-section over a broad range of neutrino 
energies, nuclear targets and Q 2 is a primary goal of the MINERvA experiment. We present 
preliminary comparisons of data and simulation in a sample rich in v^p — >> ji+n events from 
approximately one eighth of the total V events collected by MINERvA to date. We discuss future 
CN ■ plans for quasi-elastic analyses in MINERvA. 
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^ The MINERvA experiment was designed, in part, to study charged-current quasi- 

elastic (CCQE) neutrino (v^n —> /i~p) and anti-neutrino (v^p —> /i + ft) scattering of 
neutrinos with energy 1 to 10 or more GeV over a broad range of Q 2 . MINERvA 
is currently running in the "low energy" (LE) NuMI beam designed for the MINOS 
experiment in which it plans to integrate 4 x 10 20 protons on target (POT) in neutrino 
mode and 2 x 10 20 POT in anti-neutrino mode. The anti-neutrino exposure is a mixture 
of data taken with a partially constructed detector and the full MINERvA detector. After 
a shutdown beginning in 2012, the NuMI beamline will operate in the "medium energy" 
(ME) tune. Figure 1 shows the predicted neutrino flux in both the LE and ME beams. 

The active elements of the MINERvA detector [1] are polystyrene scintillator bars 
formed into detector planes. A continuous series of these planes forms the "tracker" 
region of the MINERvA detector. The analysis reported in this paper used events exclu- 
sively from the tracker; however, the MINERvA detector also contains a region where 
the planes of scintillator are interspersed with passive nuclear targets. These will allow 
study of A-dependence in the CCQE rates at different energy or Q 2 . Figure 1 also shows 
the expected rate of CCQE interactions in the scintillator tracker and in the passive nu- 
clear targets. Of course, not all of these events will be accepted in the analysis, but 
the combination of the significant mass of the MINERvA targets and the high intensity 
statistics available in the NuMI beam are very promising for this measurement. 

The planned results for CCQE in MINERvA include measurements of the CCQE total 

cross-section, <7, and da jdQ 2 on the scintillator target. At low to moderate Q 2 , Q 2 ~ 
1 GeV 2 , the primary goal is comparison with results from K2K [2, 3], MiniBoone [4, 
5], SciBooNE [6] and NOMAD [7]. In particular, there are open questions about the 
apparent cross-section enhancement at low energy and moderate Q 2 and suppression at 



1 http://minerva.fnal.gov 
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FIGURE 1. (left) The Neutrino Flux in the NuMI Low Energy (LE) and Medium Energy (ME) beams. 
Current MINERvA data is from the LE configuration, and future running concurrent with the NOvA 
experiment will use a spectrum more similar to the ME beam, (right) Fiducial CCQE Interactions in the 
Low Energy (LE) beam. The fiducial mass is also given in metric tons for each target. 



low <2 in the existing data. MINERvA will repeat these measurements on its multiple 
nuclear targets to look for changes with increasing A and nuclear density. At higher 
<2 2 , MINERvA will be able to make the first measurements of the axial form factor 
and compare the high Q 2 behavior with vector form factors. MINERvA also plans to 
compare v e CCQE with the same process using the ~ 1.5% of its events from v e 
contamination in the beam. Finally, MINERvA will search for the analogous neutral 
current process vp —> vp which is experimentally more challenging due to backgrounds 
from non-leptonic processes. 

The study described in this paper is a preliminary comparison of muon kinematics 
in CCQE candidates between data and simulation with a sample of 0.4 x 10 20 POT 
in the LE anti-neutrino beam. These data were taken during the construction of the 
MINERvA detector when only approximately one half of the detector had been installed. 
The scintillator tracker fiducial mass for this sample is 3.0 of the 6.4 metric ton fiducial 
mass available in the full detector. This sample is 12% of the total exposure measured as 
the product of POT and fiducial mass. The reference simulation is GENIE 2.6.2 [8]. The 
CCQE cross-section at nucleon level comes from the derivation of Llewellyn-Smith [9] 
with vector form factors from the BBBA2005 parametrization [10], the pseudo-scalar 
form factor from PCAC, and the axial form factor in the dipole form with ma — 
0.99 GeV/c 2 . The nuclear effects when scattering from carbon come from the Fermi gas 
model of Bodek and Ritchie [11], and Pauli blocking is implemented by a requirement 
that the outgoing nucleon has momentum of above the Fermi momentum of 221 MeV/c. 
Final state interactions of outgoing hadrons from the hard scattering process is simulated 
by INTRANUKE [8]. 



SELECTION OF CCQE CANDIDATES 



There are many possible approaches to the selection of CCQE candidates in the MIN- 
ERvA detector. We expect, eventually, to have at least three strategies for selection, each 
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FIGURE 2. (left) A candidate event in the search for v^p —> {i+n. The track originates in the scintillator 
region and is momentum analyzed in the MINOS near detector. The 30 MeV energy deposit away from 
the track is consistent with the direction of the recoil neutron, (right) For simulated events, the recoil 
energy away from the vertex vs. Q 2 derived from the muon energy and angle. The region selected in this 
space keeps very high efficiency but lets in significant background at high Q 2 . 



of which will select different CCQE-rich final states and admit different backgrounds in 
the analysis. First, we can perform a topological track based analysis with a one-track 
(muon) low Q 2 component and a two-track (muon and nucleon) component at higher Q 2 
where reconstruction of the nucleon, at least when it is a proton, is feasible. This style 
of analysis can be most easily compared to results using similar techniques from Sci- 
BooNE [6] and NOMAD [7]. Second, we can perform an inclusive analysis that selects 
only a muon and removes background based on a veto on Michel electrons from the 
chain 7i + —> /i + V jU — >> ^ + v^V jU V jU which is most similar to the analysis from the Mini- 
BooNE [4, 5] experiment. Finally, one can identify tracks as in the first analysis but 
remove events with additional visible recoil energy in the detector which is far from the 
vertex. The latter stipulation is important for keeping events with low energy nucleons 
in the final state from the remaining nucleus after the CCQE interaction. The analysis 
presented here as our first preliminary work towards this reaction is in the final of these 
three styles. As we progress in this analysis, we plan to add the other techniques. 

This analysis first selects events with a /i + originating from the scintillator tracker 
whose momentum and charge are measured in the MINOS near detector [1]. From the 
/i + reconstructed momentum and angle with respect to the neutrino beam direction we 
can reconstruct the neutrino energy and Q 2 , 

E Y tc _ m n + 2E V ( m P ~ gg) ~ ( m P ~ E b) 2 ~ m l 
2 {m p -E B -E^+ cos M ^E 2 -m 2 ^ 

and(£ 2 ) rec = 2Ef c fa -cosO^E 2 -m 2 ) -m\, 

where rn n , m p and are the neutron, proton and muon masses, E^ and 6^ are the muon 
energy and angle, and Eb is the target proton binding energy, 34 MeV in our estimation 
of the neutrino kinematics. 
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FIGURE 3. Distributions of the position at which the muon in v^p — >> ji+n originates, (left and middle) 
The transverse directions, x and y within the hexagonal fiducial volume, (right) The direction along the 
beam, including regions outside the fiducial volume selected by the analysis such as the more dense 
electromagnetic calorimeter (ECAL) region downstream of the fiducial volume and the very dense hadron 
calorimeter (HCAL) farther downstream. 



A quasi-elastic candidate event passing this selection is shown in Fig. 2. The only 
other visible energy in the event is the single bar with a 30 MeV deposit consistent with 
the expected direction of recoil based on the direction and energy of the muon. The 
expected recoil is a neutron with kinetic energy of 150 MeV. Typically only a small 
fraction of the neutron energy is observed calorimetrically in the MINERvA detector. 

For this first analysis, we chose a very conservative cut on the recoil energy in order 
to keep the efficiency of the analysis high. This approach, however, lets in significant 
backgrounds, particularly at high Q 2 as shown in Fig. 2. For events originating within 
the fiducial volume, the asymptotic efficiency at high energies is w 65%. The efficiency 
is highest at low Q 2 and falls slowly to « 30% at high Q 2 not because of the low recoil 
selection, but because of the decrease in muon acceptance to enter the MINOS near 
detector. A hit-level GEANT4-based simulation of the detector is used to evaluate these 
efficiencies and correct the neutrino interaction simulation for detector efficiencies and 
resolutions. We gain confidence in this simulation by comparison of distributions of 
candidate events in the data with the simulation in quantities that are largely independent 
of the neutrino flux and the interactions. For example, Fig. 3 shows that the distribution 
of candidate vertices in the transverse and longitudinal dimensions of the detector agrees 
with a relatively normalized simulation. 

RESULTS 

We present the results as comparisons to the simulation, including background events 
generated as described above. The simulation also includes our prediction of the untuned 
flux from a beamline simulation which has large uncertainties that we expect to reduce 
in the future [12]. We evaluate these uncertainties in the a priori neutrino flux to be 
lowest, about 7%, at the peak energy of the flux shown in Fig. 1 and significantly larger 
above this "focusing peak" energy with an asymptotic uncertainty at high energy of 
16%. These large flux uncertainties must be included in the comparisons between data 
candidates and the simulation. These uncertainties allow a significant distortion in shape 
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FIGURE 4. E v (left) and Q 2 (right) reconstructed from muon kinematics by Eqns. 1 and 2, respectively, 
for v^p —> \i+n candidates. Different reactions in the simulated event sample are shown stacked in the 
solid histogram. The flux uncertainty is shown as a band around the simulation prediction. 
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FIGURE 5. The relatively normalized ratio of data to simulated Q 2 distributions for E y < 3 GeV (left) 
and 3 < E y < 5 GeV (right). The uncertainties shown are statistical only; flux uncertainties approximately 
cancel in the analysis of the shape of this distribution, detector systematics and background uncertainties 
which are not included in the shown uncertainties may not. The change in shape expected if ma = 1 .3 GeV 
is also shown. 



of the neutrino flux with neutrino energy, E y . However, since the Q 2 spectrum as a 
function of E v is approximately independent of E v , at least far below the kinematic 

limit, Q L ~ 2m p E v , we do not expect the flux uncertainties to significantly distort the 
predicted Q 2 distribution. 

Figure 4 shows the reconstructed neutrino energy and Q 2 of the reactions for candi- 
date events. The comparison shows a disagreement outside the uncertainty band in the 
distribution of neutrino energies just above the focusing peak, 3 GeV< E v < 6 GeV. 
However, the shape of the Q 2 distribution, at least in the region where CCQE events 
make up the majority of the signal, is in reasonable agreement with our simulation. Note 
that the simulation uses an axial mass of 0.99 GeV. MiniBoonE and SciBooNE, mea- 
suring quasi-elastic neutrino and anti-neutrino scattering for 0.3 < E v < 2 GeV, prefer 
significantly higher axial masses, consistent with ma = 1.3 GeV [4, 5, 6], whereas the 
NOMAD experiment, measuring quasi-elastic neutrino scattering, 3 < E v < 100 GeV, 
prefers ma consistent with the value of 0.99 GeV assumed in our simulation. Figure 5 
shows the ratio of data to simulation for the MINERvA analysis in the low and high 



neutrino energy portions of the sample. Both appear more consistent with the low ma 
assumed rather than the high of MiniBooNE and SciBooNE; however, caution is 
warranted in drawing this conclusion because the MINERvA result only includes sys- 
tematic uncertainties from the neutrino flux, and does not yet include uncertainties on 
the detector reconstruction or the backgrounds to the quasi-elastic interactions. 

CONCLUSIONS 

The MINERvA collaboration has performed its first comparison of a quasi-elastic rich 
anti-neutrino sample with expectations for a range of assumed axial masses. The se- 
lection requires a muon from a neutrino interaction in the scintillator target and a low 
observed recoil away from the vertex in the final state to enhance elastic events. MIN- 
ERvA possesses an anti-neutrino sample approximately eight times as large as that in 
this analysis on tape, and will be able to improve the background rejection of this analy- 
sis in the future through a variety of techniques that the capabilities of the detector make 
possible. 
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